Mon. Not. R. Astron. Soc. 000, [HS] (2009) Printed 7 August 2009 (MN WI^ style file v2.2) 



Discovery of the orbital period in the supergiant fast X-ray 
transient IGR J17544-2619 



Accepted 2009 August 5. Received 2009 August 5; in original form 2009 July 16 



O ■ D. J. Clark^*, A. B. Hill^, A.J Bird^, V. A. McBride\ S. Scaringi\ A.J Dean^ 

^School of Physics and Astronomy, University of Southampton, Southampton. S017 IBJ, UK 
Laboratoire d'Astrophysique de Grenoble, UMR 5571 CNRS, Universite Joseph Fourier, BP 53, 38041 Grenoble, France 

O 



ABSTRACT 

The supergiant fast X-ray transient (SFXT) system IGR J17544-2619 has displayed 
many large outbursts in the past and is considered an archetypal example of SFXTs. A 
search of the INTEGRAL/ISGRI data, archive from MJD 52698-54354 has revealed 11 
outbursts and timing analysis of the light curve identifies a period of 4.926±0.001 days 
which we interpret as the orbital period of the system. We find that large outbursts 
occasionally occur outside of periastron and place an upper limit for the radius of the 
supergiant of <23Rq. 

Key words: Gamma-rays: observations - X-rays: binaries - X-rays: individual: IGR 
J17544-2619 
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1 INTRODUCTION 



Discovered o n the 17th Septembe r 2003 by the IBIS/IS GRI 
instrument (jUbertini et al.l 12003 Lebrun et al.l 20031 ) on- 



board the INTEGRAL satellite (IWinkler et al 
J 17544 



200 j ). 



IGR 



2619 has become the archetypal example of the 
Supersiant Fast X - ray Transient (SFXT) class of objects 
jSguera et al.ll2005l : iNegueruela et al.ll2006l ). Around twenty 
SFXTs are known to date and as the name suggests this 
class of object is characterized by rapid outbursts, lasting 
typically a few hours, with a flux 10^-10^ times the faint 
quiescent emission X-ray luminosities in the range of 10'^'^- 
10'''* ergs s~^. Through analysis of optical and NIR spec- 
tra and periodicities in the X-ray light curves, SFXTs have 
been shown to consist of a supergiant donor and a compact 
accreter. The observation of spin periods in a few SFXTs 
(IGR J16465-4507, AX J1841. 0-0535 and IGR J18483- 
0311) show that the compact object is often a neutron star. 

A number of hypotheses to explain the flaring be- 
havior of SFXTs have been proposed. A common in- 
terpretation invokes a structured wind from the super- 
giant companion, either in the form of clumping of 

a spherically symmetric outflow from th e donor star 

llii7t Zand 1005.; Walte r & Zurit a Heras 20071 : iLevder et al] 
I2OO7I : iDucci et all i2009i ) . or a wind with an eq uatorial den- 
sity p rofile inclined to the compact object orbit (jSidoli et al.l 
|2007^. In either c ase va riations in the orbital eccentricity 



I Negueruela et al] I2OO8I ) can enhance the transient nature 
of the emission from these supergiant binary systems. A 
gated accretion has also been proposed (|Bozzo et al.l 120081 ) 
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where transitions between different accretion regimes due to 
changes in the magnetospheric radius cause step-like changes 
in the luminosity. 



The initial discovery of IGR J17544-2619 on the 17th 
September 2003 was triggered by the INTEGRAL observa- 
tion of a 2 h outburst and an 8 h outburst on the same day, 
with a further outburst of 10 h detected on 8th March 2004 
demonstrating that this was a fast transient source. Subse- 
quently, more bursts fro m the source were also discovered in 
older BeppoSAX data (lin't Z and et al. 2004) and recently 
seen in XMM-Newto n IIGonzalez-Riestra ct al. 2004) . Chan 
dra 



in't Zandll2005| ) Swi. 



200£; Krimmetal 



[200 



uift 
I) 



I Krimm et"alj|2007l: ISidoU et al 
and Suzaku data fRampv et~ 



20091 ). Due to observation strategies and the relatively rare 
occurrences only a few bursts are observed a year and there- 
fore identification of any binary system parameters is chal- 
lenging. 



IGR J17544-2619 is situated ~3° from the Galactic cen- 
tre {I = 3.24°, b = -0.34°), a region shown to co ntain many 
recurrent transient sources (|Kuulkers et al.ll2007D . The loca- 
tion of IGR J17544-2619 is in a well studied region covered 
by a large amount of INTEGRAL archive data and is now 
covered by the Galactic Center Bulge Monitoring campaign. 
In this paper we report on the detection of a modulation of 
the gamma-ray fiux interpreted as the orbital motion and 
report new outbursts discovered by a search of the INTE- 
GRAL light curve. 
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2 PREVIOUS OBSERVATIONS 

The first reported de tection of IGR J175 44-2619 was on 
17th September 2003 (|Sunvaev et ahlbOOST ) when two flares 
reaching a maximum of IGOmCrab (18-25 k eV) were de- 
tecte d by IBIS/IS GRI (|Sunvaev et al.t ,2003: Grebenev et alj 
I2OO3I ). A further 5 bursts were seen when looking back at 
the older BeppoSAX-WFC data (|in't Zand et al.ll2004l . see 
Table [TJ and another from 6 d before the INTEGRAL obsex - 
vation in XMM-Newton data dConzalcz-Ri estra et al.ll2004l ). 
In total 8 bursts have been previou sly re port edly in the IN- 
TEGRAL data jWaher fc Zurita HeradfeoOTt) . ahhough the 
exact times have not been published for all the bursts. 

Observations using XMM-Newton and Chandra pro- 
vided precise positional errors of 4" and 0.6" respectively 
supportir ig the optical/NI R counterpart 2MASS J17542527- 
2619526 (|Rodriguej[200l) over the RO S AT canAiAaXe IRXS 
J175428.3-2620 (jWiinandi^ '200 ^) Using com bined optical 
and X-ray observations Pcllizza et al.l (|2006l ) classified the 
counterpart of IGR J17544-2619 to be a 09I b with a mass 
(rf25-28M0. Using optical and NIR data, iPeUizza et"al] 
l|2006t ) estimated a distance of 2. 1 -4.2 k pc based on extinc- 
tion measurements. iRahoui et all (|2008l ) used the data from 
iPellizza et al.l (|2006l ) to fit the SED and obtain a radius 
range for the donor star of 12.7Rq < 7?* < 26.6-Rq. 

An analysi s of Chandra data during quiescence by 
lin't Zandl l|2005l ') produced a power law fit with F = 5.9±1.2, 
much softer than that expected for a black hole, suggesting 
that the compact object was a neutron star. If the accretor 
is a black hole then a radio flux of ~ 10-35 mjy is expected. 
However, an upper limit of 7mJy has been placed on the 
radio emission of IGR J17544-2619 adding t o the evidence 
that the compact object is a neutron star (|Pcllizza ct al.' 
120061) . In the analysis of Suzaku observations, Jlampy ct al. 
l|2009h found that the absorption showed varied rapidly on 
timescales of minutes by a factor of ~2 during flaring, inter- 
preted as a wind clump passing though the line of sight. 
However, longer integrations showed stable val ues of Ng 
compatible with the stable absorption seen by lin't Zandl 
l|2005D . but this may simply be the result of averaging. 

The initial bursts found by INTEGRAL led to a pre- 
dictio n of a 165±3 day period in the source l|Walter et al.l 
l200d ). This period was supported by the few bursts found 
in BeppoSAX and Swift data. Long period system s can be 
difficult to explain l|Brandt fc Podsiadlowskilll995l ). A long 
period requires a large eccentricity to explain any substan- 
tial accretion, but this would decay quickly and require flne 
tuning not to disrupt the binary. Therefore, it would be un- 
likely to see too many of these systems. No evidence for 
spin or orbital periods has been observed in XMM-Newton 
or Chandra data. 



3 DATA ANALYSIS 

Archival IBIS data from MJD 52698 to MJD 54354, 
giving an on source exposure of ~ 8 Msec, were pro- 
cessed with_Jlie_J7yT£GRAi Off-line Science Analysis 
(OSA iGoldwurm et al.l l2003h version 7.0. A light curve 
for IGR J17544— 2619 was generated on science window 
timescales (~ 2000 s) over the 18-60 keV energy range. 
The light curve was then searched for bursts. The shortest 




Figure 1. Overlaying a 185 day 40% duty cycle period on the 
18—60 keV ISGRI science window light curve clearly shows that 
this period is a result of the INTEGRAL observing strategy. 



timescale region of the light curve that produces the high- 
est detection significance greater than 4(t is found. This is 
the brightest burst in the light curve. This section of the 
light curve is then removed from the total light curve and 
the next highest significant section is searched for. In to- 
tal 11 bursts have been found in the light curve (Table [1} 
which can be added to the list of bursts found in the litera- 
ture. The bursts range from short fiares of 0.5 h to periods 
of fiaring activity up to 35 hours. The time between bursts 
ranges from less than 9 days to 197 days. By folding the 
burst times on different periods, an estimation of the duty 
cycle required over the orbit to explain all the bursts can be 
made. There is a clear minimization of ~40% duty cycle at a 
period of 185 days in our data. However, this period can be 
shown to be a result of the observing strategy for the source 
(Fig. [1]). INTEGRAL is restricted in observing the Galactic 
center due to visibility constraints including the 40° solar 
aspect angle due to the fixed solar arrays and so there is an 
~6 month period to the observations of IGR J17544-2619. 
Overlaying this period and duty cycle on the ISGRI light 
curve, it is clear that this period is a result of the observing 
strategy, and with only a few points of refere nce this could 
be mistaken for a 165 day period suggested by I Walter et al.l 
(2006). No other significant periods were found by the initial 
folding analysis. 



4 PERIODICITY ANALYSIS 

The science window light curve (containing 6443 indepen- 
dent data points) was searched for signs of periodicity usin 
the Lomb-S cargle periodogram method (Lomb 1976; Scargl( 
1 19821 . Il989l ). The Lomb-Scargle periodogram is shown in 
Figure [2] with a clear peak evident at a frequency of 0.203 
days~^, corresponding to a period of 4.926 days. The signif- 
icance of the peak was confirmed by performing a random- 
ization test l|Hill et al.l [2005i) ; the fiux points of the light 
curve are randomly re-ordered and a new Lomb-Scargle pe- 
riodogram generated, the distribution of the power of the 
largest peak is an indicator of the significance of detection. 
From 200,000 simulations we estimate the 99.9%, 99.99% 
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Table 1. Bursts found in literature and in this pape r. T he peak flux is the 1 8-60 keV maximum average science window flux, assuming 
an Nh = 3.8 x 10^^ and power law index a = 2.25 l|Gonzalcz-Ric stra et al.| [2004y ★ represents bursts found in our light curve. (1) not 
seen in INTEGRAL light curve due to observation gap. ^ 0.5-10 keV. 



Date 


Orbital 


Duration 


Sig 


Peak Flux 


Instrument 


Reference 




MJD 


Phase 


Hours 


a 


lO^^crgcm^^s^^ 








50320.6 


0.39 


3.3 






SAX/WFC 


in't Zand et al. (2004) 




51229.7 


0.94 


0.2 






SAX /WFC 


in't Zand et al. (2004) 




51248.8 


0.81 


0.4 






SAX/WFC 


in't Zand et al. (2004) 




51807 5 


0.23 


8.3 






SAX/WFC 


in't Zand et al. (2004) 




51825.1 


0.80 


1.0 






SAX/WFC 


in't Zand et al. (2004) 
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0.75 


~2.0 
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53062 59 


0.01 


1.12 


5.76 


28.6±6 


INTEGRAL /ISGRI 






5S072 S5 


0.99 


5.34 


30.13 


119. 7±6 


INTEGRAL /ISGRI 


Grebenev et al. (2004) 




53189 


0.67 


0.8 




230. Ot 




in't Zand (2005) 


(1) 


53269.88 


0.09 


4.22 


10.41 


39. ±11 


INTEGRAL/ISGRI 


Seuera et al. (2006) 


* 


53441.23 


0.88 


0.56 


10.59 


62.6±6 


INTEGRAL/ISGRI 


Seuera et al. (2006) 


* 


53481.26 


0.00 


1.88 


12.20 


42.1±5 


INTEGRAL/ISGRI 




* 


53628.67 


0.93 


5.06 


6.13 


44.1±9 


INTEGRAL/ISGRI 




* 


53656.92 


0.66 


34.69 


12.99 


73.5±6 


INTEGRAL/ISGRI 




* 


53806.87 


0.10 


3.75 


16.91 


70.4±7 


INTEGRAL/ISGRI 




* 


53987.52 


0.77 


2.06 


8.32 


34.5±5 


INTEGRAL/ISGRI 




★ 


53998.38 


0.98 


4.59 


12.61 


102.4±17 


INTEGRAL/ISGRI 






54364.3 


0.26 








INTEGRAL/ISGRI 


Kuulkers et al. (2007) 




54412.1 


0.96 








Swift 


Krimm et al. (2007) 




54556.9 


0.35 


0.13 






Swift, Suzaku 


Sidoli et al. (2009). Ramov et al. (2009) 




54570.4 


0.10 








INTEGRAL/ISGRI 


Kuulkers et al. (2007) 




54708.1 


0.05 








INTEGRAL/ISGRI 


Kuulkers et al. (2007) 




54713 


0.05 








Swift 


^Sidoli et al. (2009) 




54905 


0.02 


~0.55 






Swift 


Krimm et al. (2009) 





and the 99.999% significance levels which are shown in Fig- 
ure [S] note that the maximum power achieved in any of the 
randomized hght curves was <22. 

The error on the identified period was estimated using 
a Monte-Carlo simulation. Each fiux measurement was ad- 
justed using Gaussian statistics within its individual error 
estimate to generate a simulated light curve of the source. 
The corresponding periodogram was produced and the lo- 
cation of the maximum peak near the frequency of 0.203 
days~^ was recorded. From 200,000 simulations we estimate 
the period and its error to be 4.926±0.001 (lu) days. 

Figure |3] shows the phase folded light curve using a pe- 
riod of 4.926 days and using an ephemeris of MJD 52702.9 
to centre the fold at phase 0.0. Since we would expect any 
enhancement to the emission to coincide with the neutron 
stars closest approach to the donor, we consider phase 0.0 to 
be periastron for the rest of the paper. Between phases 0.25 
and 0.85 the folded light curve is consistent with zero which 
could be interpreted as an eclipse. However, this eclipse 
would last 3d, 60% of the orbit and so is unlikely. To 
see the effect of the outbursts on the phase folded curve a 
new fold was performed excluding all of the outbursts de- 
tected (Table [!}. This second phase folded light curve is 
also shown in Figure [S] and shows no substantial difference 
in shape, demonstrating that there is an underlying peri- 
odicity in the flux aside from the bursts when then entire 
light curve is considered. The public light curves for IGR 
J1755-2619 from the RXTE- ASM and Swift-BAT missions 
were also analyzed using the Lomb-Scargle periodogram to 
see if there was any indication of the 4.926 day periodicity. 



Neither data set showed an indication of a periodic signal. 
Folding the RXTE-ASM and Swift-BAT light curves on the 
INTEGRAL period similarly showed no indication of a mod- 
ulation. 



5 OUTBURST RECURRENCE 

Aside from the bursts that have been reported in Table [T] 
there is the possibility of lower level emission occurring at 
each periastron pass that does not show up as a burst during 
our search. In order to test the recurrence of emission at 
periastron we carry out a recurrence analysis j Bird et"aLl 
1200 9). For this test we calculate the light curve significance 
for a region around the orbital period and at the midpoint 
of the period. This will then show the recurrence of bursts 
at periastron and if we are missing any low level emission 
seen at every periastron. 

In the case of IGR J 17544-26 19 the bursts do not all 
occur at periastron, reducing the sensitivity of this test. 
Since this orbit is short (~5 days) we take a region of 1 
day either side of periastron (Fig. [Jjl . Comparing the distri- 
butions of significances for apastron and periastron using a 
Kolmogorov-Smirnov test we find them to be different at the 
99% level. Out of the 79 periastrons seen only 11 produce 
a significance over 3a. The difference between the distribu- 
tions at periastron and apastron can be explained by the 
bursts seen and not an underlying longer timescale varia- 
tion as seen in SAX J1818. 6-1703 or the Suzaku data for 
IGR J17544-2619. If we repeat this with the known bursts 
removed the KS test finds that the two distributions are not 
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Figure 2. The Lomb-Scargle periodogram of the IBIS/ISGRI I8- 
60 keV light curve. The 99.9%, 99.99% and 99.999% significance 
levels derived from randomization tests are over-plotted. 



Figure 4. Recurrence analysis of IGR J17544— 2619 light curve. 
A comparison of the 2-day significance at periastron (solid line) 
and the 2-day significance and apastron (dotted line) shows that 
the only difference is duo to the bursts. 



0.5 
Phase 



Figure 3. The IBIS/ISGRI 18-60 koV Ught curve of IGR 
J17544— 2619 folded on 4.926 days and using an ephemeris of 
52702.9 days. The narrow line is for the whole data set, while the 
bold line is the result when all of the outbursts are removed. 



Figure 5. Folding the bursts on 4.926 days shows that burst 
may occur at any time during the orbit, but are more probable 
at phase 1.0, likely to be periastron. 



significantly different demonstrating the underlying emis- 
sion, which can be seen with all the phase data added to- 
gether in the Lomb-Scargle analysis, cannot be seen when 
considering individual periastrons. 

6 DISCUSSION 

This is the third SFTX with a period of a few days. 
The other two having very si milar properties to IGR 
J17544-2619. IGR J16418-4532 jCorbet et al.1 |2006| ) and 
IGR J16479-4514 jjain et al.ll2009l ) have periods of 3.75 d 
and 3.32 d respectively. This contradicts the view that 
SFXTs have different orbital geomet ries to the classical 
HMXBs (|Waher fc Zurita Her"a3l2007l ') and suggests that a 
difference in the stellar wind may lead to the difference be- 
tween SFXTs and HMXBs. 

Using the ephemeris calculated by the timing analysis, a 
histogram of burst times with respect to orbital phase can be 
produced (Figure [5]). This shows while bursts may occur at 
any time during the orbital period they are most likely in the 



last 20% of the orbit, suggesting that this phase is the point 
of closest approach of the orbit between the two bodies. 
Since we have a mass for the supergiant of 25-28M0 and 
assuming a mass of 1.4 Mq for the compact object, which is 
probably a neutron star, Keplers 3"^"^ law gives a semi-major 
axis of 36-38 Rq for the 4.9d period. T he range of possible 
radii for an Q9Ib s tar are 12.7-26.6R0 l|PeUizza et al.ll2007l : 
iRahoui et al]|2008h . 

Figure [6] shows how the Lagrange point between the 
two stars chang es with orbit for different eccentricities 
(|Paczviiskilll97ll ). We do not see any persistent or regular 
emission every orbit and so there should be no Roche-lobe 
overflow in the system. This limits the size of the star to 
<23Rq with no eccentricity in the orbit. If the minimum 
possible radius of the star is used then this limits the eccen- 
tricity of the orbit to less t han <0.4. 

In SAX J1818. 6-1703, IZurita Heras et al. I (|2009l ) pro- 
posed that the wind and compact object velocities will limit 
the region over which the NS can accrete to less than a few 
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Figure 6. The Lagrange point radius between tlie neutron star 
and supergiant as a function of phase for different eccentricities of 
orbit (e = 0, 0.1, 0.2, 0.3, 0.4 solid lines). If the Lagrange point is 
less than the supergiant radius (dashed lines are the upper/low er 
limits) then Roche-lobe overflow is expected l|Pac2vnskilll97ll) . 



stellar radii. With the short period detected in IGR J 17544- 
2619 the neutron star wil l stay relatively close to the donor 
star. IZurita Heras et al. I (j2009l ) describe how the accretion 
rate will vary strongly with distance from the supergiant, 
which may be the cause of the underlying modulation seen 
in the folded hght curve (Fig. [3]). With bursts seen at all 
phases this suggests the neutron star is within this accre- 
tion envelope and the variation could be due to the change 
in wind/NS velocity. This variation could also be due to 
eccentricity effecting the probability of seeing a burst. 

The 10'^^ ergs s ~^ quiescent emission seen by Suzaku 
l|Rampv et al.l |2009| ) is very low and not easily explained 
by steady, spherical ly symmetric Bondi-Hoyle accretion 
l|Bondi fc Hovlel 19441 ) which gives a mass loss rate for the O- 
type star of M ~ 10"^° Mq yr"^ based on a NS at ~2Rv,. 
The increase in luminosity to 10^^ ergs s~^ during a burst 
would require a 4 magnitude increase in the accretion rate, 
bringing the mass loss rate up to the more common 10~^ 
Mq yr~^ usually seen from a supergiant donor. This change 
cannot be explained simply by the change between perias- 
tron and apastron. One interpretation would be that the 
wind is clearly not homogeneous and isotropic. If we sug- 
gest that the bursts are caused by a clumpy wind or a sim- 
ilar outflow, the probability of the neutron star interacting 
with a clump that keep s a constant size as it moves out will 
be proportional to 1/r^. lDucci et al] (|2009l ) showed that the 
size of a clump changes as it moves away from the donor 
star, Kci (x (r^(l — where /3 is a constant in the 

range ~0.5-1.5 dependent on the supergiant and for the rest 
of this paper assumed to be 1.0. This will result in the prob- 
ability of interacting with a clump being proportional to 
{r'^^^ — r^^^)/r'^ resulting in a change in the probability even 
with a slight eccentricity in the orbit (Figure (Tj . If we take 
the normalized probabilities under these curves between 0.8 
to 0.2 at periastron and 0.3 to 0.7 at apastron we can com- 
pare these probabilities calculated from the distribution of 
phase folded bursts. We observe 18 bursts at periastron and 
2 at apastron out of a total of 23 bursts. This results in 
relative probabilities of 0.8 and 0.1 for the phases of 0.8 to 



Figure 7. If the probability of a stellar wind clump hitting the 
neutron star is proportional to a (r^^^ _j-i/3)/f.2 j^w then even a 
slight eccentricity (e = 0, 0.1, 0.2, 0.3, 0.4 shown) can cause this 
probability to change greatly around the orbit. 

0.2 and 0.3 to 0.7 respectively, although the low number of 
observed bursts, especially at apastron, will skew the proba- 
bilities in favour of more bursts being seen at periastron. The 
low statistics in our sample for this letter currently prohibit 
any realistic attempt at testing this model. However, with a 
larger sample of bursts from a large number of systems the 
models for a clumpy structured wind could be tested. 



7 CONCLUSIONS 

We have shown that IGR J17544-2619 has an orbital pe- 
riod of '^4.9d resulting in a semi-major axis of 36-38 Rq. 
As there is no evidence of Roche-lobe overflow, this puts an 
upper limit on the size of the donor star of ~23Rq with 
no eccentricity in the orbit. However, some eccentricity is 
required to explain the variation of emission around the or- 
bit and so the supergiant radius must be closer to its lower 
limit of 12.7Rq. With a larger sample of bursts and with 
the techniques explained here, it may be possible to better 
constrain the eccentricity and therefore the supergiant ra- 
dius, and lead to a much better understanding of the SFXT 
systems seen by INTEGRAL. 
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